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LONG-TERM  GOALS 

The  long-term  goal  of  this  effort  is  to  improve  the  Navy  community  ocean  circulation  model 
ROMS/TOMS  by  incorporating  astronomical  tidal  forcing  and  the  latest  developments  in  turbulent 
mixing. 

OBJECTIVES 

The  principal  objective  of  this  research  is  to  improve  subgrid-scale  parameterization  in  Navy 
community  and  operational  ocean  circulation  models.  This  is  to  be  accomplished  by  assessing  and 
refining  turbulent  mixing  parameterization  as  well  as  including  comprehensive  direct  astronomical 
tidal  forcing  of  importance  to  many  semi-enclosed  marginal  seas. 

APPROACH 

This  project  complements  extremely  well,  the  AESOP  program,  the  ONR  DRI  on  subgrid-scale 
parameterization  and  skill  assessment  of  numerical  ocean  models  as  well  as  the  new  Characterization 
and  Modeling  of  Archipelago  Strait  Dynamics  (CMASD)  DRI.  Using  the  Adriatic  Sea  ROMS/TOMS 
as  the  test  bed,  we  will  incorporate  direct  astronomical  forcing  of  the  1 1  major  tides  in  the  global 
ocean:  semidiurnal  M2,  S2.  N2,  K2;  diurnal  Kl,  01,  PI  and  Ql;  long  period  Mf,  Mm  and  Ssa.  The  co¬ 
oscillating  baro tropic  tides  will  be  prescribed  from  LHK’s  tidal  model  of  the  Mediterranean  Sea  (see 
ocean.colorado.edu/~kantha).  Note  that  many  global  ocean  models  such  as  the  ones  resulting  from 
NASA  initiatives  do  not  perform  well  in  some  marginal  seas  and  it  is  essential  that  a  regional  model  be 
used.  Note  also  that  compound  tides  such  as  M4  will  be  generated  by  the  nonlinear  model  itself  and  is 
an  indirect  result  of  the  principal  astronomical  tidal  forcing. 

The  latest  Kantha  and  Clayson  (2004)  turbulent  mixing  model  based  on  second  moment  closure  will 
be  incorporated  into  ROMS/TOMS.  This  model  includes  the  effect  of  surface  waves.  No  other 
turbulence  model  does  at  present.  The  inclusion  of  surface  wave  effects  such  as  wave  breaking  and 
Stokes  production  should  greatly  improve  the  simulation  of  drifter  trajectories  in  the  Adriatic. 

A  very  high-resolution  version  of  ROMS/TOMS  will  be  set  up  for  the  Venice  Lagoon  and  will  be 
driven  at  the  boundary  by  the  Adriatic  Sea  ROMS/TOMS.  This  will  facilitate  testing  the  skill  of  this 
model  in  shallow  lagoon  systems  and  is  highly  pertinent  to  the  current  focus  of  the  ONR. 
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WORK  COMPLETED 


Observational  data  to  compare  with  turbulence  models  are  scarce.  Microstructure  measurements  have 
not  become  a  routine  staple  of  oceanographic  measurements  as  CTD  casts  have  been  for  decades.  In 
collaboration  with  Dr.  Sandro  Camiel  of  ISMAR,  Italy,  who  has  a  related  NICOP  grant  from  ONR,  we 
have  participated  in  NATO  Undersea  Research  Center/Naval  Research  Laboratory  (NURC/NRL) 
DART  06A  and  06B  cruises  in  March  and  August  of  this  year,  and  collected  turbulence  data  using  a 
microstructure  profiler  (e.g.  Figure  1).  The  five  hundred  and  twenty  microstructure  profiles  collected 
greatly  enhance  the  turbulence  database  in  the  Adriatic,  which  consists  of  approximately  73  profiles 
collected  by  Peters  and  Orlic  (2005)  and  Peters  et  al.  (2006)  in  the  Northern  Adriatic.  Along  with 
other  complementary  oceanographic  data  collected  during  the  cruises,  these  micro  structure 
measurements  provide  a  nice  dataset  to  assess  turbulent  mixing  parameterization  in  ROMS/TOMS. 

See  Camiel  et  al.  (2006)  for  more  details. 

This  project  started  in  February  22nd  of  this  year  and  already  we  have  made  considerable  progress  in 
the  modeling  arena.  With  the  help  of  our  Italian  and  American  colleagues,  the  ROMS/TOMS  model 
has  been  successfully  ported  to  our  Sun  workstation,  configured  for  the  Adriatic  Sea  and  numerous 
runs  have  been  made.  The  version  we  are  running  incorporates  General  Ocean  Turbulence  Model 
(GOTM)  developed  by  EU  Turbulence  researchers,  thanks  to  Dr.  John  Warner  of  USGS.  This 
provides  a  good  framework  for  making  modifications  to  the  turbulent  mixing  parameterization  in 
ROMS.  For  example,  we  have  included  non-local  mixing  and  Stokes  production  effects  into  the  LMD 
mixed  layer  model  component  of  ROMS  and  have  investigated  the  resulting  differences. 

The  upcoming  tasks  for  this  year  include  assessment  and  refinement  of  turbulence  parameterization  in 
ROMS/TOMS  by  comparison  with  DART  data.  The  latest  Kantha  and  Clayson  (2004)  turbulent 
mixing  model  based  on  second  moment  closure  will  be  incorporated  into  ROMS/TOMS. 

We  will  work  on  incorporating  astronomical  tidal  forcing  into  ROMS/TOMS  in  the  coming  year. 

Using  the  Adriatic  Sea  ROMS/TOMS  as  the  test  bed,  we  will  incorporate  direct  astronomical  forcing 
of  the  1 1  major  tides  in  the  global  ocean:  semidiurnal  M2,  S2.  N2,  K2;  diurnal  Kl,  01,  PI  and  Ql; 
long  period  Mf,  Mm  and  Ssa.  The  co-oscillating  barotropic  tides  will  be  prescribed  from  LHK’s  tidal 
model  of  the  Mediterranean  Sea  (see  ocean.colorado.edu/~kantha). 

RESULTS 

A  destabilizing  buoyancy  flux  at  the  ocean  surface  leads  to  convective  mixing  in  the  water  column. 
Under  pure  convection,  the  TKE  dissipation  rate  £  must  simply  scale  as  the  surface  buoyancy  flux 
JbO.  It  has  been  the  practice  hitherto,  following  Shay  and  Gregg  (1984,  1986),  Lombardo  and  Gregg 

£  ~  cJ 

(1989)  and  Brainerd  and  Gregg  (1993  a&b)  to  assume  that  the  dissipation  rate  /,()  is  constant  in 
the  entire  mixed  layer  under  pure  convection  (e.g.  Peters  et  al.  1988,  2006).  The  value  of  the  constant 
is  taken  as  -0.58  following  Lombardo  and  Gregg  (1989).  However,  Camiel  et  al.  (2006)  show  that  a 
more  reasonable  value  for  c  to  be  0.39.  Therefore,  in  the  convective  mixed  layer 


F  —  J 

=  0.39,7, 


b  0 


=  0 


2  =  0 

0.1  D<z<  0.9 D 
z>D. 


(1) 


On  the  other  hand,  when  the  turbulence  in  the  mixed  layer  is  mechanically  driven,  by  the  wind  stress, 
the  law  of  the  wall  demands  that  the  dissipation  rate  near  the  surface  follow  the  relationship 


£  — 

*  ^  where  K  is  the  von  Karman  constant,  u*  is  the  friction  velocity  and  z  is  the  distance  from 
the  surface.  This  similarity  relationship  should  hold  in  the  upper  few  meters  near  the  surface  if  we 
ignore  the  wave  effects  on  8  scaling.  The  falling  microstructure  probe  did  not  allow  us  to  make 
measurements  in  the  upper  2-3  m,  where  the  influence  of  surface  waves  on  the  TKE  dissipation  rate  is 
most  prominent.  Carniel  et  al.  (2006)  show  that  in  the  wind  stress-driven  mixed  layer, 


ss  =  uI/{kz)  0  <  z  <  0.3  D 

=  3.33 uH(kD)  0.3 D<z<D 


=  0  z>D. 
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Figure  1:  Profiles  of  temperature  (°C)  and  salinity  (psu),  density  (kg  m-3)  and  buoyancy  frequency 
(s-1),  and  the  Thorpe  scale  (top  panels),  TKE  dissipation  rate  (W  kg-1),  temperature  variance 
dissipation  rate  (K2  s-1),  eddy  diffusivity  K  (m2s-l),  and  the  heat  diffusivity  Kh  (m2s-l)  (bottom 
panels)  as  measured  during  OP  B90-4  at  Station  B90  in  the  Gulf  of  Manfredonia,  under  moderate 
winds  and  nocturnal  cooling.  A  total  of  32  casts  were  made  over  2.5  hr  centered  around  the 
midnight  of  August  24th/25th.  The  green  line  denotes  the  corresponding  average  value.  The  slow 
change  in  salinity  during  the  OP  is  due  to  the  ship  drifting  over  a  patch  of  brackish  water.  Thick 

redline  in  panel  1  shows  mean  salinity. 


Below  the  mixed  layer  and  in  the  interior  of  the  water  column,  mixing  is  episodic  and  internal  wave 
field-driven.  The  relevant  length  scale  is  the  Ozmidov  length  scale 


(3) 


If  we  further  assume  that  the  Thorpe  scale  LT  ( Thorpe  1977)  is  proportional  to  the  Ozmidov  scale  LO 
(e.g.,  Dillon,  1982;  Stansfield  et  ai,  2001),  the  dissipation  rate  can  be  taken  to  be 

£,  =  0.03LpV3 


where  the  proportionality  constant  has  been  determined  by  the  best  fit  to  values  appropriate  to  the 
observed  background  dissipation  rate  deep  in  the  water  column  (depth  ~  60-80m). 
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Figure  2.  Observed  dissipation  rates  observed  at  Station  B90  compared  with 
theoretical  scaling  for  OPs  B90-2,  B90-3  and  B90-4:  thick  black  line  (Eq.  5),  black 
dotted  line  (Eq.  6),  and  green  line  (observational  mean). 

When  the  turbulence  is  generated  by  both  the  momentum  flux  and  a  destabilizing  buoyancy  flux,  the 

TKE  dissipation  rate  in  the  mixed  layer  can  be  taken  to  be  the  sum  of  the  rates  due  to  shear-driven 
and  buoyancy-driven  turbulence.  Therefore 


z<D 
z>  D 


+  SS 


=  £; 


(5) 

Figure  2  shows  the  TKE  dissipation  rate  profiles  plotted  along  with  the  profile  indicated  by  Eq.  (5)  for 
OPs  B90-2,  B90-3  and  B90-4.  The  conventional  scaling  ( Lombardo  and  Gregg  1989,  Brainerd  and 
Gregg  1993  a&b,  Stips  et  al.  2002) 


sc  =  0.58 Jb0;  ss  =  1 .76i/*  / (atz) 

E  Sc  +  Ss 


(6) 


is  also  shown.  It  can  be  seen  that  Eq.  (5)  is  a  better  depiction  of  the  dissipation  rates  in  the  deep  than 
the  traditional  formulation  (Eq.  6),  which  has  no  validity  below  the  upper  mixed  layer  and  hence 
should  not  be  applied  except  in  the  mixed  layer.  In  the  mixed  layer  itself,  the  difference  between  the 
two  formulations  is  small,  although  Eq.  (5)  is  better  justified  from  first  principles.  The  disagreement 
between  the  theoretical  fonnulations  and  the  observed  values  is  undoubtedly  due  to  inaccuracies  in 
inferring  JbO  and  u*  from  bulk  formulae. 


IMPACT/APPLICATIONS 


Accurate  depiction  of  many  quantities  of  interest  to  worldwide  naval  operations,  such  as  the  upper 
layer  temperature  and  currents,  requires  accurate  simulation  of  turbulent  mixing  in  the  water  column 
and  accurate  tidal  forcing.  Operationally,  this  contributes  to  better  counter  mine  warfare  capabilities 
through  better  and  more  accurate  tracking  of  drifting  objects  such  as  floating  mines.  Other  drifting 
materials  such  as  spilled  oil  are  also  better  tracked  and  counter  measures  made  more  effective.  Other 
applications  include  search  and  rescue. 

RELATED  PROJECTS 

1.  Subgrid-scale  Parameterization  in  3-D  Ocean  Models:  The  Role  of  Turbulent  Mixing  (PI  -  Dr. 
Sandro  Carniel  of  ISMAR,  Venice,  Italy)  -  NICOP. 

2.  Improving  the  Skill  of  Ocean  Mixed  Layer  Models  (PI  -  L.  Kantha)  -  N00014-05-1-0759.  Ended 
June  2006. 
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